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■ Abstract: We derive effective actions for SU(2) Polyakov loops using inverse Monte 
Carlo techniques. In a first approach, we determine the effective couplings by requiring 

', that the effective ensemble reproduces the single-site distribution of the Polyakov loops. 
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The latter is flat below the critical temperature implying that the (untraced) Polyakov 
loop is distributed uniformly over its target space, the SU (2) group manifold. This allows 
for an analytic determination of the Binder cumulant and the distribution of the mean- 



field, which turns out to be approximately Gaussian. In a second approach, we employ 
novel lattice Schwinger-Dyson equations which reflect the SU(2) x SU (2) invariance of the 
functional Haar measure. Expanding the effective action in terms of SU(2) group characters 
makes the numerics sufficiently stable so that we are able to extract a total number of 14 
couplings. The resulting action is short-ranged and reproduces the Yang-Mills correlators 



very well. 
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1. Introduction 

The deconfmement phase transition in pure Yang-Mills theory [1, 2] is controlled by the 
dynamics of the Polyakov loop variable Above a critical temperature T c , the singlet 
part L x = tr develops a nonvanishing vacuum expectation value (VEV). In this high- 

temperature phase one expects to find a plasma of liberated gluons (and, in QCD, also 
quarks). The VEV of L x thus represents an order parameter associated with spontaneous 
symmetry breaking. The symmetry in question is a global Zjy symmetry, Zat being the 
center of the gauge group SU(N). While the Yang-Mills action is center symmetric, L x , 
although gauge invariant, transforms nontrivially, L x — > zL x , z £ Z^r. Combining renor- 
malization group ideas and dimensional reduction, Svetitsky and Yaffe have conjectured 
that finite-temperature SU(N) Yang-Mills theory in d dimensions is in the universality 
class of a Zjy spin model in dimension d— 1 [3, 4]. For some recent and rather sophisticated 
confirmations of the statement on the lattice the reader is referred to [5, 6, 7, 8]. 

The universality argument implies that effective field theory methods may be put 
to use. It should make sense to map the microscopic theory, here Yang-Mills, onto a 
macroscopic one, described by an effective action with Zat symmetry. For gauge group 
SU(2), for instance, one can try to coarse-grain the gauge fields all the way down to Z2 
Ising spins [9, 10, 11]. An intermediate procedure is to establish an effective action for the 
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Polyakov loop variable itself [12, 13, 14]. This may be achieved analytically using strong- 
coupling or, equivalently, high-temperature expansions [13, 15, 16]. Doing so for SU(2), 
one obtains a local effective action depending on all characters XjiVx) [15, 16]. The index 
j € INI/2 labels the irreducible representations of SU(2). In this most elementary case, Xj 
can be expressed in terms of powers of L x (the character of the fundamental representation, 
j = 1/2). For larger gauge groups, however, more and more characters/representations 
become relevant. This fact has recently been employed for model building, regarding the 
untraced holonomy ^ [17] or, equivalently, its eigenvalues [18] as the fundamental degrees 
of freedom. We parametrize the (lattice) effective action as follows, 

S C S = ^aS a , (1-1) 
a 

with center-symmetric operators S a and effective couplings A a to be determined. As stated 
above, for SU(2) it is sufficient to work with only the traced Polyakov loop, L x . The 
effective action will then have the form [14], 

S e fi[L x ] = V[L X ] + L x K X yL y + L x L y K X y UV L u L v + . . . . (1-2) 

x xy xyuv 

The kernels depend on the couplings A a and the temperature. By construction, 

the 7.2 center symmetry (L x — > —L x ) is manifest. Note that the representation (1.2) 
is rather general and leaves room for a plethora of operators, the compact continuous 
variable L x G [—1,1] being dimensionless. Later on, it will therefore be crucial to choose 
an appropriate subset of all possible operators in order to capture the essential physics. 
In this respect it turns out useful to follow [17] and view the effective action (1.1) as 
being embedded into a 'sigma model' depending on S c g[L] = 5 c ff[^Pa:]- This yields an 
additional global SU{2) symmetry, 

Vx^gWxg- 1 , gesu(2), (1.3) 

which is a remnant of the underlying SU (2) gauge invariance. The Haar measure V^ x has 
an even larger symmetry, namely SU{2) x SU(2), corresponding to the transformation law 

Vx^gVxh, g,heSU(2). (1.4) 

The invariance of the measure leads to novel Schwinger-Dyson identities which will be an 
important ingredient in our derivation of the effective couplings A a inherent in (1.1). 

The paper is organized as follows. In Section 2 we derive exact (lattice) Schwinger- 
Dyson equations from the invariance of the Haar measure T>y} x . We proceed by analysing 
the single-site distribution of the Polyakov loop variable L x in Section 3. This yields 
a semianalytic method to determine all couplings A a apart from the one of the hopping 
term, Ao- The latter is obtained in Section 4 using the Schwinger-Dyson equations which 
are also employed to check the resulting effective action. In Section 5, we determine the 
effective potential in the symmetric phase from the single-site distribution. Finally, in 
Section 6, we perform an extensive numerical analysis to improve the effective action by 
including a maximum number of 14 operators. Some technicalities concerning the analysis 
of histograms are relegated to an appendix. 
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2. Haar measure and Schwinger-Dyson identities 



The Polyakov loop variable on the lattice is given by a holonomy or parallel transport 
connecting the (periodic) boundaries in temporal direction, 

N t 

Vx = Hu t , x . , (2.1) 
t=i 

where the C/'s are the standard link variables on a lattice of size N t x N 3 . The effective 
action for the Polyakov loops is obtained by inserting unity into the Yang-Mills partition 
function, such that (the trace of) (2.1) is imposed as a constraint, 

Z YM = JvU exp(-S w [U]) 

, N t 

= / VUV^5(ii^ x -tv\[U t , x . fi ) exp(-S w [U]) 
J t=i 

= j Pqjex P (-5 efr [q}]), (2.2) 

with DU and D*p the appropriate Haar measures (see below) and Sw the standard Wilson 
action. Of course, the integration over link variables U in the last step cannot be performed 
exactly. For this reason one has to resort to effective actions as given by (1.1) and (1.2), 
for instance [3, 4, 14]. Using inverse Monte-Carlo (IMC) techniques, it should be possible 
to determine a reasonable effective action from Yang-Mills configurations. 

The main ingredient for this procedure are the Schwinger-Dyson equations associ- 
ated with the symmetry of the measure 2?*p under (1.4). To derive those we choose the 
parametrization, 

<P^P°l + 2T a P^P£^ , (2.3) 

which is in SU(2), tytcVx = 1, if the components P x define a three-sphere S 3 according 
to 

P« = (P°) 2 + P^ = 1. (2.4) 

We mention in passing that the points x where the Polyakov loop is given by center 
elements, ty x = ±1, correspond to the positions of monopoles in the Polyakov gauge 
[19, 20, 21], a particular realization of 't Hooft's Abelian projections [22]. 

In terms of the coordinates (2.3), the traced Polyakov loop becomes L x = while 
the functional Haar measure can be written as 

V<# = \[d i P x 8{P£P£-\) . (2.5) 

X 

Obviously, this is invariant under rotations R G SO (A) generated by the angular momenta 

LZ = -i(p>JL-PzJL). (2.6) 
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These can be split up into 'electric' and 'magnetic' components (or 'boosts' and 3d 'rota- 
tions'), 

iK% , (2.7) 
ie abc L c x , L a x = ^e abc L b x c . (2.8) 

Summarizing, the SO (4) generators L x v rotate the four-vector P£, while the SO (3) gen- 
erators L% rotate the three-vector P x . The self- and anti-selfdual combinations, 

MZ = ±(L%-KZ), (2.9) 

K^l&l + K), (2-10) 

generate left and right multiplication, respectively, 

Vx^gVx, ¥x^¥xh, g,heSU(2). (2.11) 

Global SU (2) (gauge) transformations of the Polyakov loop as given by (1.3) are generated 
by L x b (or L x ) which do not differentiate with respect to the trace P x and thus leave any 
functional of P x = L x invariant. Typical such invariants are 

P°, P£P« = 1-P£P°,... . (2.12) 

The Schwinger-Dyson equations that follow from the SO (4) invariance of the Haar measure 
(2.5) are given by 

J I£{F[% exp(-S eff [<p])} = , (2.13) 

where F[ty] is an arbitrary functional of ty x - As the effective action depends on *X) X solely 
through the SU{2) invariant P°, S cS [^\ = 5 eff [P°], only the generators L% a = K% lead to 
nontrivial relations which can be written as 



x dP% 

■jab _ pa 

x ~ x dPZ 



pa g 

dp® 
d 



-P b 



dP£ 



(KF[%] - FlWKZSesW]) = , (2.14) 
using the expectation value notation, 

(O) = Z- 1 J Om exp(-S cS m) . (2.15) 

Because K x transforms like a vector under gauge rotations, (2.14) in general will not be 
gauge invariant. However, we are still free to choose the functional F[ty] at our will. If we 
pick 

F^m = P^G[P°], (2.16) 
with an arbitrary functional G[P°], we have the action of K x , 

K a x F b y = 5 ab 5 xy P®G - P a x P b y G' x , (2.17) 
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where we have denoted G' x = dG/dP x . Plugging this into (2.14), setting x = y and taking 
the trace one finds the gauge invariant Schwinger-Dyson equations, 

<3P°G - P«P«{G' X - GS' eS>x )) = . (2.18) 

The same result is obtained using = K X H[P°] instead of (2.16) and identifying 

H' x = —G x . Let us rewrite (2.18) as a functional integral, 



/ 



[3P°G - P a x P a x {G' x - GS' eS;X )] exp(-SW) = , (2.19) 
and parametrize ty x according to 

= exp ir a 9 x = 1 cos 9 X + ir a n a x sin 9 X , n% = P%, '(P b x P b x ) l/2 . (2.20) 
Then, the traced Polyakov loop is L x = cosO x while the Haar measure (2.5) becomes 

UV = flam 2 (2.21) 

X 

As the functional integral (2.19) only depends on invariants we can integrate over the 
directions n (yielding an irrelevant volume factor) so that we are left with an integral 
involving only the reduced Haar measure, 

VL = Y[ d6 x sin 2 6 X = Yl d(cos 6 X ) sin 9 X = ]J dL x ^l-L 2 x = ]J DL X , (2.22) 

XX X X 

namely, 

' Yl d9 y sin 2 9 y [3 cos 9 X G - sin 2 9 X (G' X - GS' cS J] exp(-S cS ) = . (2.23) 
v 

A more compact form for these relations is achieved in terms of total derivatives, 
= J J] d9 y sin 2 9 y J d9 x — {sin 3 9 X G exp(-S cff )} 

y+x " J - 

= J \{ de y sin2 e v J d ( cos e *) » g g n ( sin3 ^ G eM-Scs)} ■ (2.24) 

Note that the sin 3 9 term ensures the absence of surface terms. With (2.24) we have found 
the Schwinger-Dyson relations of the reduced theory involving only the invariant L = cos 9. 
We do not have a simple geometrical explanation for the invariance of the reduced Haar 
measure VL leading to (2.24). The 50(4) symmetry of the measure Vty, however, is very 
natural. 

In terms of the Polyakov loop L x , (2.23) is the expectation value 

(3L X G - (1 - L 2 X ){G' X - GS')) = . (2.25) 



/: 
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Comparing with (2.18) we notice that it does not matter whether the expectation value is 
taken with the full or reduced Haar measure as long as G = G[L]. If we insert the ansatz 
(1.1), the Schwinger-Dyson equations (2.25) become a linear system for the couplings A„, 

£<(1 - L x ) GS' a:X ) \ a = ((1 - L x ) G' x ) - 3 (L X G) . (2.26) 

a 

To solve this unambiguously we need as many independent operators G as there are cou- 
plings A a . A particularly natural procedure, which also turns out to be rather stable 
numerically, is to choose G = S' b . Any of these operators contains an odd number of L x s 
so that the minimal set of Schwinger-Dyson equations relates only nontrivial expectation 
values, 

E^ 1 " L D S 'b,y S 'a, x ) X a = ((1 - Ll) S' b \ yx ) - 3 (L x S' bty ) . (2.27) 

a 

At this stage, keeping x and y fixed, the problem of determining the couplings A a is well 
posed mathematically. Numerically, of course, it is better to use all the information one 
can get, for instance by scanning through all possible distances x = \x — y\, x < N s /2. 
The resulting overdetermined system is then solved by least-square methods. Another 
possibility is to add new equations to (2.27) by choosing further appropriate monomials 
or polynomials in L x for the operator G. This philosophy will be extensively adopted in 
Section 6. Before that, however, we will try to proceed in a (semi-) analytical fashion. 

3. Single-site distributions of Polyakov loops 
3.1 Definitions 

From the effective action of Polyakov loops S e g [L] one can derive new probability densities 
by integrating over (part of) the loop variables L. Of course, this amounts to some kind of 
course-graining so that via the new densities one will only have access to gross properties 
of the effective action. Nevertheless, these densities, if chosen properly, exactly repro- 
duce certain expectation values calculated within the full effective ensemble. Consider, for 
instance, the local moments, 

i p = {Ll) = Z- 1 j l[DL y L x exp(-S eS [L}) , (3.1) 
J y 

where, as usual, the partition function Z is the integral over exp(— S^s )■ Splitting off the 
L^-integration, (3.1) can be rewritten as 

l p = (Ll) = J 1 DL x LP xPw [L x ] = (L x ) w , (3.2) 

with the probability density pw obtained via integrating over all L y / L x , 

PW [L X ] = Z- 1 J ll DL y exp (S cS [L y ]) = Z' 1 exp (-W[L X \) . (3.3) 
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Due to translational invariance, pw (like £ p ) does not depend on the site x. Thus, DLpyy[L] 
is the probability to find the value of the Polyakov loop in the interval [L,L + dL]. The 
Z2-symmetry of the effective action implies that the power p in (3.1) and (3.2) has to be 
even, p = 2q, at least for finite volume (no spontaneous symmetry breaking). Therefore, 
knowing pw gives access to all local moments li q and (by taking the logarithm) to all local 
cumulants C2 q as well. A particularly important quantity is the Binder cumulant [23, 24], 
defined as the quotient 

&4-| = |-3, (3.4) 

which measures the deviation from a Gaussian distribution. This will be analysed in some 
detail later on. 

From the definition (3.3) it is obvious that pw is blind against spatial correlations 
of Polyakov loops. In other words, one cannot calculate two-point functions like G xy = 
(L x L y ). In principle, this can be remedied by a slight generalization of (3.3). To this end 
we define a new probability density depending on L x and L y , 

pw^L x ,L y \ = Z~ x J H J DL z exp(-S , cff [L]) = Z- 1 exp(-^ 2 [L a; ,L y ]) . (3.5) 

z^x, y 

Then, one can calculate the following two-point correlators, 

(LILD = j DL x DL y LlLlp W2 [L x ,L y ] . (3.6) 

Obviously, pw and pw 2 are related according to 

Pw[L x ] = J DL y p W2 [L x ,Ly} . (3.7) 

If there were no correlations, one would have factorization, pw 2 [Lx, L y ] = pw[L x ]p\y[L y ]- 
3.2 Determination of single site distributions 

At first glance, there seems to be not much of a gain by introducing densities like the 
single-site distribution pw- Note, however, that pw[X] is much simpler than our original 
density ps = Z~ x exp(— S e g) which depends on variables rather than just one. In 
addition, pw can be obtained rather easily from our Monte Carlo data. The results are 
fairly smooth histograms which are displayed in Figure 1 (for details see App. A). The most 
important observation, however, is the finding that pw is flat below T c , that is, one has 
an equipartition for L x . Apparently, this is a remnant of the SO(4) symmetry discussed in 
Section 2. Taking the (negative) logarithm of pw we obtain the single-site potential W[L] 
shown in Figure 2. 

We are thus led to employ the following ansatz for the potential W from (3.3), distin- 
guishing between temperatures below (— ) and above (+) the critical value, T c , 



W-[L} = const, (3.8) 

2k 



W + [L] = const' + £fr^ • ( 3 - 9 ) 
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Figure 1: Single-site density for temperatures above (x,*) and below T c (+). For T <T C 

(/3 < C ~ 2.299), the density is flat, p w = 2/tt. Input: 200 to 400 configurations, N s = 20. 

Demanding (1) = 1 these imply for the density pw, 

PwlL] = exp(-I47_)/Z_ = 2/vr , (3.10) 
p+[L] = eM-W+[L])/Z + . (3.11) 

Things are particularly straightforward below T c , so let us discuss this case first. The 
result (3.10) shows that, after normalization, the single-site distribution of Polyakov loops 
below T c is known exactly. Furthermore, it is simple enough so that the associated (local) 
moments can be determined analytically, 

* - - § /_; - -rife*. (, 12) 

The generating function for these moments can also be calculated explicitly, 

Z_(t) ^ (e tL ) w _ = I [ DL e tL = £ ^U* = | J x (t) , (3.13) 
^ ^ z>o ^ 

/i being the standard modified Bessel function. For the Binder cumulant (3.4) we thus 
find the result 



6 * = 7FjS- 3 =(^- 3 = - 1 - 
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Figure 2: The single-site potential W[L] shifted by the constant offset w. For T < T c ((3 < j3 c 
2.299), W[L] is flat. Input: 200 to 400 configurations, N s = 20. 



We have checked that (3.12) and (3.14) hold numerically both for the histograms pw and 
the effective Yang-Mills probability density ps- The results for the Binder cumulant are 
displayed in Figure 3. 

It may seem strange that we get a flat distribution pw below T c . However, this does 
not imply that the effective potential, which defines the distribution of the mean field L, 
becomes trivial (see Sect. 5). 

To proceed, we have to specify our ansatz for the effective action beyond (1.1) and 
(1.2). Svetitsky and Yaffe have argued [3, 4, 14] that, close to the phase transition, the 
effective interactions should be short ranged so that 5 e ff is of Ginzburg-Landau type, 

S c s = Xo^2L x L x+i + ^2^2^(L x ) 2k ^X So + X 2 S 2 + ... . (3.15) 

The high-temperature character expansions mentioned in the introduction yield additional 
hopping terms of the form L p x L q y . . . [15, 16, 25]. The relevance of these terms will be 
dicussed in Section 6. 

Let us investigate the consequences of the ansatz (3.15) for the single-site distribution. 
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Figure 3: The Binder cumulant 64 as obtained from the simulated Yang-Mills configurations (□) 
with 7V S = 20 compared to the single-site distribution pw (*)■ Below (3 C ~ 2.30, the exact result 
(3.10) for pw (i.e. W = const) has been used. Above /3 C , W has been fitted to a polynomial (see 
below) . 

Plugging the former into the definition (3.3) yields 




where, in the second line, we have introduced the field 

n a 

M x = — ^ = Y^(L x+i + L x „i) , (3.17) 

I 

representing the sum of all nearest neighbors of L x . In addition, we have defined a modified 
action S' eS which is obtained from S e g by setting L x = 0, 

S' eff [L] = S eS [L}\ Lx=0 . (3.18) 

Now, the left-hand side of (3.16) is 2Z^/tt and hence independent of L x . Thus we may 
put L x = everywhere on the right-hand side yielding the identity, 




(3.19) 
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Accordingly, e~ w ~ is the partition function associated with action S' s . We can go one 
step further and expand the exponential containing the nearest-neighbor field M x on the 
right-hand side of (3.16). This is actually a hopping-parameter expansion in Ao which, 
upon using (3.19), implies 

Here, we have defined modified expectation values associated with S' eS and Z', 

(0[L\Y = J J] DL y O[L] eM-S' cS [L])/Z> . (3.21) 

The Z2~symmetry of the effective action requires n to be even, n = 2m. Denoting 

» 2m = (Ml m )' , (3.22) 



we finally have 



(2m)! x 1 \ ^ 2k 

m=0 y ' \k=\ 



To lowest order in L x (m = 0) this consistently reproduces the normalization (3.19), 
(1)' = 1 = e~ w ~ /Z' '. A general interpretation can be given as follows. To have equipartition 
requires a delicate balance between the hopping term (Ao) and the 'potential' terms (A2fc). 
Setting Ao = (so that the effective action leads to a product measure) implies that all 
A2fe have to vanish and vice versa: \2k = implies Ao = 0. 

To further evaluate the identity (3.23) we note that it can be viewed as a particular 
example of a linked-cluster or Mayer expansion [26, 27, 28] expressing the moments X^^m 
in terms of the cumulants 

X' 2k = (2k - 1)! A 2fc . (3.24) 

The relation between moments and cumulants can actually be solved for arbitrary m (see 
e.g. [29]), 

n=l k 1 ,...,k n =l v ' x ' 1=1 

fc 1 +... + fc n — m 

This somewhat clumsy formula yields for the first few orders 

Ao> 2 = A 2 , (3.26) 

A^ 4 = A 4 + 3A 2 2 , (3.27) 

Ao Me = A' 6 + 15 A' 2 A 4 + 15 A 2 3 , (3.28) 

Aq Ms = A' 8 + 28 A 2 Ag + 35 A 4 2 + 210 A 4 A 2 2 + 105 A' 2 4 . (3.29) 

It is quite obvious that by inverting (3.25) we can express the couplings A2fe (or cumulants 
X' 2k ) in terms of the moments fi2m- Alternatively, one may take the logarithm of (3.23) 
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and compare coefficients. In any case, the first few cumulants are 

A' 2 = A^ 2 , (3.30) 

A4 = A^(/i 4 -3^) , (3.31) 

A 6 = A^ (fiQ - 15 /x 4 ^2 + 30 nl) , (3.32) 

Ag = A(j (^ 8 -28^ 6 ^2 + 420 ^4^2 - 630^-35^4) . (3.33) 



These identities almost solve our problem of determining S e s as they express the unknown 
couplings A 2 fc in terms of Ao (unknown as yet) and the modified expectation values //2m 
from (3.22). 

Things become simple if one allows for only a finite number (say K) of couplings \2k 
in the Svetitsky-Yaffe action (3.15). Then, there is only a finite number of independent 
moments ii 2 k-> k = I, . . . K . This is quite obvious from e.g. (3.33). Setting \$ = = A' g 
determines the moment fj,$ and all higher ones in terms of 112, M4 and fj,Q. 

For K = 1, (3.23) yields the general expression 

/A \ m 

t x 2m = (2m-l)\\ i^J =(2m-l)!! A i5 l , m=l,2,.... (3.34) 

We thus have found factorization: all higher moments fi2m, m > 1 can be expressed in 
terms of the lowest one, /j, 2 = A2/AQ. Of course, this is consistent with S e g being quadratic 
in L x (vanishing of quartic and higher cumulants A^,). 

For K = 2, we have three couplings, Ao, A2 and A4. In this case, (3.23) implies the 
following generalization of (3.34), 



M2m — 




(3.35) 



which shows that all moments fi 2 m can be expressed in terms of /j, 2 and ^4. The first 
two factors in the sum count the number of ways in which one can form k pairs out of 
m elements. The term raised to power k is actually (one third of) the Binder cumulant 
associated with the moments fi2m- If it were zero we would get back at (3.34). 

Clearly, in order to determine the couplings A 2 fc one does not want to calculate the 
moments fi 2 k by performing a new and costly Monte Carlo simulation with the action S' eS , 
setting L x = at a particular site x. One expects, however, that, for large lattices, one 
will have the approximate identity 

(M 2m )' ~ (M 2m ) , m > , (3.36) 

where the latter expectation is taken in the full Yang-Mills ensemble. For our numerical 
evaluation we have tested assumption (3.36) as follows. Define the expectation values 

(M 2 x m )A = Z A 1 I H DL y Ml™ exp ( - S cS [L y ] - AL 2 X ) , (3.37) 
J y 

so that one has 

(M 2m } = (M 2m } , {Ml m )' = {M x m ) 00 . (3.38) 
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If (3.36) is to hold then (M x m }\ must be approximately independent of A. We have checked 
this by simulating the leading-order action, 

for different values of A on a lattice of size 16 3 with Ao = —0.3 (symmetric phase). The 
calculated expectation values (M 2 }\ displayed in Table 1 show that (M 2 )\ is indeed inde- 
pendent of A to an accuracy of about 0.5 %. 



A 


1 10 100 1000 10000 




1.951 1.947 1.962 1.954 1.939 1.961 



Table 1: The expectation value (M x }\ as a function of the parameter A suppressing the single-site 
variable L x . Input parameters are N s — 16, A = —0.3 (symmetric phase). 



For T > T c , we use the ansatz (3.9). This implies that formulae (3.20-3.33) still hold, 
however, with A2A; now replaced by X 2 k — «2fc- We have checked that the identification 
(3.36) also holds in the broken phase (choosing Aq = — 1, see Table 2). 



A 



(M 2 x )a 



10 



100 1000 10000 



18.79 18.87 18.83 18.78 18.80 18.78 



Table 2: The expectation value {M x )\ as a function of the parameter A suppressing the single-site 
variable L x . Input parameters are N s = 16, Ao = —1 (broken phase). 

The couplings K 2 k can be obtained by fitting W+[L] (see Figure 2) according to (3.9). 
The fit values are displayed in Tables 3 and 4. 



p 


k 2 /2 


K 4 /4 


2.40 


-0.4468 


0.0703 


2.34 


-0.2712 


0.0526 


2.32 


-0.1772 


0.0261 


2.30 


-0.0717 


0.0120 



P 


k 2 /2 


K 4 /4 


K 6 /6 


2.40 


-0.4531 


0.0901 


-0.0152 


2.34 


-0.2626 


0.0249 


0.0216 


2.32 


-0.1612 


-0.0259 


0.0408 


2.30 


-0.0666 


-0.0087 


0.0133 



Table 3: Two-parameter fit to W+[£]. 



Table 4: Three-parameter fit to 



Summarizing we note that we have good analytical and numerical control of the single- 
site distribution pw or, equivalently, the histograms displayed in Figure 1. Below T c , the 
histogram is flat, = const, above T c , W + ~ \ogp^ is a simple polynomial in L 2 with 
coefficients given in Tables 3 and 4. 
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4. Determination of the effective action 



The calculation of the couplings \2k, k > 0, in the effective action proceeds in three 
steps. First we determine the moments \i2m from the Polyakov-loop ensemble using the 
approximate identity (3.36). Second, from (3.30-3.33), we obtain the couplings \<ik = 
X' 2k /(2k — 1)!, k > 0, in terms of the moments \i2k an d Ao- Third, we determine Ao- 

The first step consists of straightforward numerics based on our Wilson ensembles 
obtained for several values of (3 near (3 C . The results for the //2m are displayed in Table 5. 






2.20 


2.25 


2.28 


2.29 


2.30 


2.32 


2.34 


2.40 




1.93 


2.086 


2.242 


2.327 


2.466 


2.946 


3.336 


4.173 


M4 


10.16 


11.55 


13.07 


13.89 


15.27 


20.16 


24.22 


33.60 


M6 


80.88 


96.06 


113.0 


121.7 


137.6 


194.1 


241.5 


357.6 


M8 


829.3 


1019 


1237 


1341 


1551 


2297 


2922 


4536 



Table 5: The moments //2m for different values of the Wilson coupling (5 (N s = 20). 
With the moments //2m at hand we find the couplings 

A 2 fc = \ 2 o k a 2 k , k > , (4.1) 

where the ct2k can be expressed in terms of the (i2k according to (3.30-3.33). The final 
step consists in the determination of Ao- To this end we make use of the Schwinger-Dyson 
relations (2.26) choosing the operators G = L^ 1 which results in 



({\-Ll)M x I%- l )\ Q +Y,{^-L 2 x )L 

k>0 



2k+2l-2 
x 



}A 2 , = (2l-l)((l-^)eV3(^). (4.2) 



-2l\ 



For T < T c , where the single-site distribution is known exactly, the right-hand side of (4.2) 
vanishes. This can either be inferred from the analytical result (3.12) or by noting that 
the term in question is a total derivative, 



(2i-i)((i-^ir 2 )-3^) 



-2l\ 



vr dL 



(1- £2)3/2 L 



21-1 







(4.3) 



Plugging (4.1) into (4.2) and dividing by Ao (assumed to be nonzero) yields a nonlinear 
equation of degree 2k— 1 in Ao- With the coefficients «2fc and all nonlocal expectation values 
(correlators) determined numerically, the coupling Ao can be obtained straightforwardly. 
As there are 2k — 1 solutions we take the one which is approximately independent of the 
number K of couplings \2k- The resulting values of all couplings (for K = 2 and K = 3) 
are displayed in Tables 6 and 7. 

With the effective couplings determined we are in the position to check our results by 
simulating the effective action. For both (3 = 2.20 and (3 = 2.40 we have produced 10000 
configurations distributed according to S e g using the couplings from Table 7. 
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p 


2.20 


2.25 


2.28 


2.29 


2.30 


2.32 


2.34 


2.40 


Ao 


-0.438 


-0.473 


-0.500 


-0.509 


-0.630 


-0.685 


-0.697 


-0.725 


A 2 /2 


0.186 


0.233 


0.280 


0.301 


0.453 


0.603 


0.675 


0.873 


A4/4 


-0.002 


-0.003 


-0.005 


-0.007 


-0.019 


-0.053 


-0.088 


-0.212 



Table 6: Numerical values for the couplings Ao and A2fc/2fc, k < K — 2. The critical Wilson 
coupling is (3 C = 2.299 (N s = 20). 



p 


2.20 


2.25 


2.28 


2.29 


2.30 


2.32 


2.34 


2.40 


Ao 


-0.438 


-0.476 


-0.507 


-0.510 


-0.628 


-0.690 


-0.705 


-0.760 


A 2 /2 


0.186 


0.237 


0.288 


0.303 


0.453 


0.621 


0.698 


0.979 


A4/4 


-0.002 


-0.003 


-0.006 


-0.007 


-0.020 


-0.057 


-0.093 


-0.256 


A 6 /6 


0.000039 


0.00011 


0.00027 


0.00037 


0.0020 


0.0106 


0.0244 


0.116 



Table 7: Numerical values for the couplings Ao and A2fc/2fc, k < K = 3. The critical Wilson 
coupling is f3 c = 2.299 (N s = 20). 



0.85 



0.75 



" 0.7 



0.65 



0.6 



0.55 



-0.6 



-0.4 



-0.2 




L 



0.2 



0.4 



YM * 
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K = 3 x 



0.6 



0.8 



Figure 4: Comparison of single-site histograms based on simulating Yang-Mills (*) vs. the effective 
action for T < T c . The curves for two and three couplings A2fc, i.e. K = 2 (+) and K = 3 (x), 
respectively, fall on top of each other. Input: (3 — 2.20, N s = 20. 
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0.55 1 1 1 1 1 1 1 1 1 1 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

L 

Figure 5: Comparison of single-site histograms based on simulating Yang-Mills (*) vs. the effective 
action with two (+) and three (x) couplings \2k for T > T c . Input: j3 = 2.40, N s = 20. 



p 




Ao A2 A4 A6 


2.20 
2.20 


input 
output 


-0.43803 0.37182 -0.00681 0.00024 
-0.43824 0.37351 -0.00621 0.00020 


2.40 
2.40 


input 
output 


-0.76000 1.9572 -1.0216 0.69761 
-0.76027 1.9605 -1.0222 0.69039 



Table 8: Comparison of couplings used as input of simulation with couplings obtained as output 
of IMC applied to the effective action. 



In Figures 4 and 5 we compare the single-site distributions obtained from the effective 
theory with those of Yang-Mills. The outcome is quite satisfactory. In particular, one 
notes that the inclusion of a L 6 -term (K = 3) still improves the matching of the histograms 
compared to the case K = 2. 

A further important check is provided by reproducing the input couplings of Table 7 
via our IMC procedure. The results displayed in Table 8 show quite convincingly that the 
method works. If we allow for additional operators in the numerics (which are not present 
in the effective action) the numbers of Table 8 remain unchanged while the couplings of 
the new operators are consistently of order 10 -5 , i.e. compatible with zero. 
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5. The constraint effective potential 



With an effective action being found, one could go on and calculate the constraint effective 
potential [30] which defines the distribution of the constant mean field, 

Z ^E L *' ^ = iV s 3 . (5.1) 

X 

In perturbation theory, the effective potential has been evaluated long ago [31, 32]. It 
describes a 'gas' of gluons at high temperature, i.e. deep in the deconfined phase. Recent 
models for the effective potential which also describe the confined phase are based on the 
eigenvalues of the Polyakov loop ty x [18] and not just their sum L x . As stated in the 
introduction, this difference becomes obsolete for 577(2). 

It thus seems of interest to investigate the effective potential on the lattice. This 
apparently requires further Monte-Carlo simulations of the effective action S e s [L] with the 
mean field L held fixed, following the approach adopted in [30, 33]. It turns out, however, 
that these additional efforts can be avoided by making use of some statistical properties of 
the single-site distribution p\y discussed in Section 3. 

The constraint effective potential V is defined in terms of the probability density of 
the mean field (5.1), 

p v [L] = Z^e-WW = Z- 1 I VL5(l-SI- 1 Y,Lx) exp(-5 cff [L]) , (5.2) 

J X 

with the normalization Zy given by the partition function 

Zy = Z v (0) = J 1 dle- nv ^ . (5.3) 

In what follows, we will try to obtain the mean-field distribution py from the single- 
site distribution p\y- We note, first of all, that, due to translational invariance, the first 
moments coincide, 

(L)y = j dLL Pv [L] = Y, Z 1 III DL V Lx e ~ ScS[L] = = W • ( 5 - 4 ) 

•> x -> y 

The higher moments, on the other hand, are different, 

(L p )w = [ilDL y LI e- s °«W = {L P) (5 . 5) 
J y 

(L*>)v = n-P Yl (L xl ...L Xp )= X ^ . (5.6) 

For the mean-field distribution we thus get generalized susceptibilities while pw yields 
expectation values of arbitrary powers of L at a single spatial site, taken in the ensemble of 
Polyakov loops extracted from Yang-Mills. This has been discussed at length in Section 3. 
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To obtain a connection between arbitrary moments we suppose that the generating 
functions associated with py and pw are related according to 



Z v (t) = (exptZ)v = ([[exp(tL x /Q)) v ~ J](exp(iL/ft)) w = [Z w (t/n)f . (5.7) 

X X 

Here, we have made the assumption that only a small fraction of the random variables 
{L x : x e ft} are statistically dependent. This is justified for large volumes and short- 
range correlations. According to the law of large numbers we expect the collective random 
variable L = J2 X to have a Gaussian distribution if the L x are randomly distributed 1 . 

Let us check to which extent this is realized. 

Below T c , Zw = Z- is exactly known from (3.13) so that 



Zvit) 



2ft i n _ W 



t -h(t/n) 



J.ZK 

= Ej^yCL 2k ), Z V (0) = 1. (5.8) 



Thus, by expanding the Bessel function (to power ft) we know all moments or susceptibil- 
ities of py • Explicitly, one finds 

<JV = ± ■ <5 ' 9) 
^ - w + 3J w^ • (5 ' 10) 

f6 5 15(0-1) 15(fl - l)(fi - 2) 

{L)v 64g + 32^ + 64S? • (5 - U) 

In the large-volume limit, ft — > oo, the leading terms yield 

^ 2 V = ^^ = (2A;-1)!!(L 2 )^, (5.12) 

an identity typical for a Gaussian distribution. As a cross check, we calculate the Binder 
cumulant associated with py. From (5.9) and (5.10) we have 

^ = ^- 3 — s- (513) 

which obviously vanishes in the infinite-volume limit in accordance with (5.12). Summing 
up the moments (5.12), we obtain the large-volume partition function 

Z v (t) ~ exp(t 2 /8ft) , (5.14) 

which turns out to be Gaussian in t. Substituting t = iu, we have 

Z v (iu) = J dl exp(-ftV[Z] + iuL) ~ exp(-u 2 /8ft) . (5.15) 

1 Note, however, that with L being a compact variable, we cannot expect a Gaussian in a strict mathe- 
matical sense. 
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0.25 



Figure 6: Gaussian fits to the distribution py[X] obtained from simulating Yang-Mills on lattices 
of size 16 3 x 7 and 7 3 x 7. The value Nt = 7 for the temporal extension corresponds to the symmetric 
phase. 

To extract the mean-field distribution py = exp(— QV)/Zy we take the Fourier transform 
with respect to u and find 



p v [L] ~ y/2n/ir exp(-20L 2 ) , (5.16) 

which is a perfect Gaussian distribution with variance 

a 2 = 1/40 = (L 2 ) v . (5.17) 

The fact that L is compact does not really matter as in the large- volume limit assumed, 
the Gaussian is sharply localized at L = 0. This is indeed seen from Figure 6 which shows 
that a Gaussian fit to the distribution of L, 

Pvfit [L] = -=L- exp(-L 2 /2<r 2 ) , (5.18) 

V ITXO 



works perfectly well. 

This is corroborated by comparing the fit values for a with the expectation values 
calculated from Yang-Mills as displayed in Table 9 for different volumes and bin sizes. 

The agreement between the fitted width and the expectation value (Z 2 ) 1 / 2 is quite 
impressive, in particular for large volumes, as expected. Due to the approximations made, 
however, we do not reproduce the absolute numbers given by (5.17). If we define 
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n x N t 


config.s/bin 


a 


(L 2 ) 1 / 2 


7 3 x 6 


120 


0.0837 


0.0773 


7 3 x 6 


80 


0.0845 


0.0773 


7 3 x 7 


120 


0.0582 


0.0549 


7 3 x 7 


80 


0.0588 


0.0549 


16 3 x 6 


250 


0.0249 


0.0252 


16 3 x 7 


150 


0.0167 


0.0164 


16 3 x 7 


250 


0.0167 


0.0164 



Table 9: Width a of the Gaussian fit (5.18) compared to the expectation value (Z 2 ) 1 / 2 calculated 
from the SU(2) Monte Carlo ensemble. The values for the temporal extension N t correspond to 
the symmetric phase. 

we get for Oi = 7 3 and Q2 = 16 3 the numerical value 7 = (16/7) 3 = 11.94 while the results 
of Table 9 yield 

7 = 11.1 ±0.4, N t = 6, (5.20) 
7 = 11.5 ±0.9, N t = 7, (5.21) 

where the error has been estimated by varying the bin sizes. Thus, at least for sufficiently 
low temperature (large Nt) we obtain the correct scaling of the width with the volume. 

6. Reproducing the two— point function 

The procedure developed so far is based on the single-site distribution of the Polyakov loop 
which is under good (semianalytic) control. By construction, the effective action obtained 
in this way reproduces the Yang-Mills distribution quite well (recall Fig.s 4 and 5). At 
this point it is natural to ask how well we are reproducing correlators of the Polyakov loop. 
After all, these are intimately related to the confining potential (T < T c ) or the Debye mass 
(T > T c ), see e.g. [14]. In Fig.s 7 and 8 we compare the Yang-Mills two-point function 
with the one obtained from the Svetitsky-Yaffe effective action (3.15) using the couplings 
from Table 8. 

The figures suggest that we are doing quite well in the symmetric phase ((3 = 2.20, 
i.e. T < T c ). In the broken phase (j3 = 2.40, i.e. T > T c ), however, there is room for 
improvement both in the exponential decay and the value (L) 2 of the plateau. To assess 
the (dis)agreement quantitatively, we fit all two-point functions according to 

G x0 = (L X L ) = a[exp(-bx) + exp ( - b(N s - x))] + c . (6.1) 

The values for the fit parameters are listed in Table 10 and corroborate the qualitative 
statements made above. 
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Figure 7: The Yang-Mills two-point function (YM) compared to the one obtained from the 
Svetitsky-Yaffe effective action with four couplings (sim). Input: (3 = 2.20, N s = 20. 



0.3^ 
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2 4 6 8 10 12 14 16 18 20 

X 

Figure 8: The Yang-Mills two-point function (YM) compared to the one obtained from the 
Svetitsky-Yaffe effective action with four couplings (sim). Input: (3 — 2.40, N s = 20. 

In order to improve the matching between the effective theory and Yang-Mills we 
obviously have to include more operators. In previous applications of IMC, this has mainly 
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a b c 


2.20 


YM 
sim 


0.2493 1.9627 0.0009 
0.1971 1.8309 0.0001 


2.40 


YM 
sim 


0.2006 2.0715 0.0752 
0.1295 1.4499 0.0802 



Table 10: Comparison of the fit parameters from (6.1) associated with Fig.s 7 and 8. 

been done for Ising systems [10, 11, 34, 35] or twodimensional nonlinear sigma models 
[36, 37]. In these cases, the set of operators is restricted as they square to unity. For the 
Polyakov loop, however, the situation is different, as arbitrary (ultralocal) powers as well as 
hopping terms associated with arbitrary powers are allowed, i.e. terms like L^L^L^ — 
It turns out the the IMC procedure tends to get destabilized upon including more and 
more monomials in L x . As a result, the values for the couplings depend rather strongly 
on the number of operators present and of equations used in the overdetermined linear 
system. In addition, the determinants of the matrices to be inverted may become as small 
as 10~ 40 . We thus had to work with symbolic programs like Maple, setting the number of 
digits to 60 or even more. Nevertheless, the instabilities prevailed. Inspired by the results 
from the high-temperature expansion on the lattice [15, 16], we have tried to overcome 
these problems by changing our operator basis from monomials in L to characters. Being 
orthogonal class functions, these seem to be the natural candidates for an economic set of 
operators. At this point it should be noted that for an effective action with a finite number 
of terms different choices of bases are not equivalent. 

As stated in the introduction, for SU{2) the characters can be expressed as polynomials 
in the traced Polyakov loop, L = tr*p/2 = cos 9, according to 

^'^^'Ehks:;) 1 ^-^' ,-0,1,1.... . <a,) 

This formula allows to reobtain the L-representation from the characters. The first few 
relations are 

Xi/ 2 = 2L , Xl = 4L 2 - 1 , xs/2 = 8L 3 - 4L , .... (6.3) 

These are sufficient to obtain monomials up to terms of order L^L y . To streamline notation 
it is useful to define a basic link variable associated with lattice points x and y and SU (2) 
'color spin' j, 

X j]Xy = Xj(L x )Xj(Ly), (6.4) 



- 22 - 



• — • « // * I — • f/—^ • — • I m 

-0.11150 -0.02003 -0.00477 0.00257 0.00368 0.00191 -0.00052 
-0.15908 -0.06020 -0.00614 0.00649 0.00535 0.00547 0.00003 



Z K 



K 



« • 



0.00090 -0.00085 0.00070 -0.00004 0.00021 -0.00833 0.00008 
0.00096 -0.00053 0.00052 -0.00055 0.00001 0.04305 -0.00061 



Table 11: Effective operators and couplings for (3 = 2.20 (upper entries) and (3 = 2.40 (lower 
entries), N s = 20. 



which we represent graphically as 



X 



1/2; xy 



x. 



l;xy 
3/2; xy 



AL x L y , 
16L 2 x L 2 y -4L 2 x 



4L y + l 



6444 



y 



32L x Ly + \QL x Ly , 



(6.5) 
(6.6) 
(6.7) 



A link with n 'internal' lines thus corresponds to the representation labelled by j = n/2. 
These links are the basic building blocks of our basis of effective operators. The leading 
order of the high-temperature expansion [15, 16] is then given by the nearest-neighbor 
expression, 

~~ (6.8) 



<SlO = ^jXj;x,x+i 1 

x,ij 



with Xj a known function of the temporal Wilson coupling (3t and extension Nt that de- 
creases rapidly with 'color spin' j. If we rewrite the basic link (6.4) as Xj- x ^ x+r , we have 
two parameters controlling our basis, the representation label j and the effective range 
('link length') r = \r\. Several test runs of the IMC routines have confirmed good conver- 
gence in j so that we will restrict ourselves to the lowest representations. The maximum 
range we allow for is the plaquette diagonal, i.e. r < To further restrict the number of 
operators, we limit ourselves to a maximum number of four links of type (6.4) that can be 
drawn within a single plaquette. A typical term, for instance, is thus given by 



K 



= X\/2- x - x +i X\/2- x,x+j Xi/2- t x,x+i+j Xl/2; x+i,x+j 



(6.9) 



Altogether we have 14 operators corresponding to 18 monomials in L. They are displayed in 
Table 11 together with the couplings associated with them. Several comments are in order 
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Figure 9: The Yang-Mills two-point function (YM) compared to the one obtained from the 
character action with 14 couplings (sim). Input: [3 = 2.20, N s = 20. 



at this point. By allowing for all possible distances x = 0, 1, . . . , 10 in the Schwinger-Dyson 
equations (2.27), we obtain a maximum number of 140 equations for our 14 operators. The 
values of the couplings remain fairly stable if we vary the number of equations used in the 
IMC least-square routine (changes being approximately 1% for the relevant couplings). 
The x 2 P er degree of freedom is 5 • 10 -5 for the maximum number of 140 equations. 

For the operators (6.5 - 6.7) we find rapid decrease with spin j. Similarly, if we 
increase the number of links within the elementary plaquette, the associated couplings 
tend to decrease. The leading order hopping term, • — • (r = 1) dominates by one order 
of magnitude compared to the terms with r = \pl. This already indicates that the effective 
interactions are short-ranged in accordance with the Svetitsky-Yaffe conjecture. 

If we enumerate the couplings by gi, . . . from left to right, we may express the new 
effective action as 

14 



>cff 



(6.10) 



a=l 



Note that, according to (6.5 - 6.7), the old LO coupling Ao is given by a (rapidly convergent) 
series in j, 



Aq = 4<?i + 16 <?i4 + terms with j > 3/2 



-0.445 for = 2.20 
-0.646 for (3 = 2.40 



(6.11) 



- 24 - 



0.3 



0.25 



0.2 

o 

0.15 



0.1 



2 4 6 8 10 12 14 16 18 20 

X 

Figure 10: The Yang-Mills two-point function (YM) compared to the one obtained from the 
character action with 14 couplings (sim). Input: [3 = 2.40, N s = 20. 








a b c 


2.20 


YM 
sim 


0.2493 1.9627 0.0009 
0.2509 1.9837 0.0001 


2.40 


YM 
sim 


0.2006 2.0715 0.0752 
0.2051 1.9257 0.0758 



Table 12: Comparison of the fit parameters from (6.1) associated with Fig.s 9 and 10. 

These numerical values for Ao agree reasonably well with those of Table 8, where 
only four operators had been used. The benchmark test to be performed, however, is the 
calculation of the two-point function G x q using the new effective couplings g a . Fig.s 9 and 
10 show that we have indeed improved the matching between Yang-Mills and the effective 
action. 

This is quantitatively confirmed by repeating the fits of (6.1) and Table 10. The 
results displayed in Table 12 convincingly show the improvement in the effective action, in 
particular for the broken phase. 
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7. Summary and discussion 

In this paper we have derived effective actions describing the dynamics of the (traced) 
Polyakov loop variable L x = tr ^3^/2, and hence of the deconfinement phase transition. It 
has turned out useful, however, to regard the effective action as being derived from a more 
general theory depending on the untraced Polyakov loop ty x [17]. This theory is a nonlinear 
sigma model with target space SU (2) = S 3 and hence the symmetry SU (2) x SU (2) 
corresponding to left and right multiplication of tyx by group elements. Although the 
effective actions in L clearly do not have this symmetry, it is nevertheless inherited by 
the functional Haar measure which implies novel Schwinger-Dyson equations for Polyakov 
loop correlators. In addition, it seems that a remnant of this symmetry shows up in the 
single-site distribution p\y of L x which is flat below T c meaning that tyx is distributed 
uniformly over the group manifold. Obviously, it would be desirable to really prove this 
equipartition for which we have found convincing numerical evidence. As the single-site 
distribution of L x is exactly known in the confinement phase, we can give exact predictions 
for all moments (L 2k ) and for the Binder cumulant, 64 = —1. Above T c , we have fitted 
the log-distribution W ~ logpw by polynomials so that also in this case we have good 
quantitative control of the distribution. 

It turns out that W[L] and a Ginzburg-Landau (or Svetitsky-Yaffe) effective action 
<Seff[L] are related in a manner that is simple enough to proceed by analytic means. Assum- 
ing that expectations taken in the effective action are unchanged if L x is changed at a single 
site (another relation valid numerically but still subject to a proof) we have been able to 
express the effective couplings \2k {k / 0) of 5 e fr[L] in terms of the parameters of W. The 
remaining coupling Ao is then determined by means of the Schwinger-Dyson equations. 
The single-site distributions resulting from the effective theory S c s [L] agree very well with 
those obtained directly from Yang-Mills. Furthermore, the Svetitsky-Yaffe effective action 
perfectly fulfills the Schwinger-Dyson equations based on the 50(4) invariance of the Haar 
measure. 

For the symmetric phase (T < T c ) we have also determined the (constraint) effective 
potential from the single-site distribution pw assuming that the interactions are sufficiently 
short-ranged such that the law of large numbers may be invoked. As expected we obtain 
a Gaussian distribution for the mean field L if the volume is large and the temperature 
small enough. 

By definition, one cannot calculate correlations from single-site distributions. Vice 
versa, the matching of these distributions does not imply that the correlation functions 
match as well. A direct comparison shows that the two-point functions of the Yang-Mills 
and Svetitsky-Yaffe ensembles differ somewhat, in particular in the broken phase. To im- 
prove the matching we have changed our operator basis from monomials in L x to characters, 
which are orthogonal polynomials in L x . Technically, this results in a numerically rather 
stable inverse Monte Carlo procedure, even if the number of operators is large. We have 
obtained the effective couplings for a total number of 14 operators. The resulting effective 
theory has short-range interactions and reproduces the Yang-Mills two-point function in 
both phases very well. 
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Further research will be devoted to the following issues. The predictions of the effective 
actions for the dynamics of the phase transition should be investigated in detail. This 
includes an analysis of the effective potential(s) near and beyond the transition point as 
well as calculations of critical exponents. The latter will yield a check whether the effective 
action S e ff[L] is indeed in the universality class of the Z2-Ising-model. In addition, it 
should be possible to generalize the methods developed in this paper to higher SU(N) 
gauge groups. Work in these directions is under way. 
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A. Histograms and bins 

Given a probability density pw[L] one defines the associated (cumulative) distribution 
function 

r L 

P W [L] = J dLWl-L' 2 p w [L'} . (A.l) 

Density pw and distribution P\y are related to our histograms as follows. We have a total 
number N of 'events' or 'measurements' saying that a Polyakov loop at site x belonging to 
an arbitrary configuration takes its value in some prescribed interval ('bin'). Accordingly, 
A is a fairly large number, 

N = Ng x A config = 20 3 x 400 = 3.2 x 10 6 . (A.2) 

The number of bins (labeled by integers i) is denoted by I, the number of events in bin 
i by Cj. This number represents the height of the ith column in the histogram counting 
the absolute numbers of events with values in [Lj_i,Lj] . The relative counting rate is 
obtained by normalization, 

d = Ci/N = P w [U] - P w [Li_i] = Pw [Li] \Jl-L* AL, , (A.3) 

where ALj = L; L — Lj_i and L; L G [Xj-i, Li] chosen appropriately. The situation is depicted 
in Figure 11. 

Good statistics is achieved if the counting rate c% is approximately constant because 
then all bins will be equally 'populated'. This can be achieved by suitably choosing the bin 
sizes ALj which, however, is somewhat tricky because of the nontrivial measure in (A.l). 
If we ignore this for the moment and choose an equidistant partition, 

ALi = AL = 2/1 , (A.4) 

the total count in bin i becomes 

Ci=™p w [Li} y/l-L?. (A.5) 
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Figure 11: General histogram for distribution function P^[L]. 



This yields rather bad statistics near the boundaries L = ±1, in particular for T > T c , 
due to the suppression by the measure. For instance, choosing [3 = 2.4, / = 100, i.e. 
AL = 1/50, one typically finds C\ ~ 14000 data points in the first bin (near L = — 1), 
while the population of the bins near L = is larger by a factor of five, C50 — 73000. The 
suppression by the geometry thus 'wins' against the density which is peaked near L = ±1. 
In the quantity of interest, the probability density, 

pwfc] = -£= J- , (A.6) 

one divides by the measure factor which tends to zero near L = ±1. This yields the 
peaks near L = ±1 but at the same time further enhances the statistical error close to 
the boundaries. For T < T c , this is not much of a problem as we have equipartition, 
Pw[L] = const = = 2/ it, and the density is known anyhow. For T > T c , however, (A.6) 
implies that the bulk of the density is located where the statistical error is largest. On 
the other hand, the behavior of pw in this regime determines the higher order couplings 
K 2 k- The lesson to be learned is that the partition should be modified such as to correctly 
incorporate the effect of the measure. To this end, we demand that the counting rate be 
constant, Cj = c, for T < T c , hence, from (A. 3), 

c = p w yJl-L? ALi = 1/1 . (A.7) 

Thus, in order to properly take into account the measure, the bin size ALi has to be chosen 

such that 

\fl-Lf ALi = const = c/p^ = —— . (A.8) 

This can be achieved by going over to continuum notation, 

c/Pw = r dL ^- L2 = p w[Li\ - Pw[Li-i] , (A.9) 

JLj_i 
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and solving this recursion for Lj numerically with Pw(L) given by 

P W [L} = ^ L Vl - L 2 + arcsin(L) . (A.10) 

Alternatively, one may produce an ordered list of all data points for L, and partition this 
list in such a way that all bins contain the same number C of 'events'. The sampling points 
Li are then given by the smallest (or largest, depending on the counting convention) value 
of L in bin i. 

For T > T c , the density pw is then given by 

Pw[U] = Ij ^C i . (A.ll) 

This has been displayed in Figure 1. Obviously, measure effects are now absent and the 
difference between Cj and C represents the deviation from equipartition. 
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